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Abstract 

We present the results of a Monte Carlo simulation of the RP^ 
model in three dimensions with negative coupling. We observe a sec- 
ond order phase transition between the disordered phase and an anti- 
ferromagnetic, unfrustrated, ordered one. We measure, with a Finite 
Size Scaling analysis, the thermal exponent, obtaining v = 0.784(8). 
We have found two magnetic-type relevant operators whose related rj 
exponents are 0.038(2) and 1.338(8) respectively. 
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The theory of Critical Phenomena offers a common framework to study 
problems in Condensed Matter Physics (CMP) and in High Energy Physics 
(HEP). In both areas, the concepts of Spontaneous Symmetry Breaking 
(SSB) and of Universality allow to relate problems in principle very different. 

The usual Heisenberg model, associated with the standard ferromagnetic 
Non Linear a Model (NLcxM), has a SSB pattern of type 5*0 (3) /SO (2). 
With the introduction of nontrivial AntiFerromagnetic (AF) interactions 
the SSB pattern normally changes completely and, usually, frustration is 
generated. In particular, a SSB pattern shared by several AF models is 
S0{?>) X S0{2) / S0{2). For instance, some frustrated quantum AF Heisen- 
berg models or the Helimagnets and Canted spin systems [|, Q are 
examples of this behavior. Frustrated quantum spin models are specially 
interesting because of their possible relation with High Temperature Super- 
conductivity Q]. 

As a general consequence of the Weinberg Theorem the low energy physics 
of a system is completely determined by its SSB pattern, the effective La- 
grangian for the system being the corresponding NLaM. In this framework 
a study has been carried out for 50(3) x S0{2) / S0{2) in perturbation the- 
ory , where the main conclusion reached is that the only possible nontrivial 
critical point in three dimensions is that of 0(4). In spite of that, we have 
found in a nonperturbative lattice formulation of the same model, a critical 
point with exponents clearly different from those of 0(4). 

From the HEP point of view, it is of great interest to understand whether 
AF interactions can generate new Universality Classes. One could even hope 
that nontrivial antiferromagnetism would be the ingredient needed in order 
to nonperturbatively formulate interacting theories in four dimensions. 

In a previous work we found that, on a three dimensional AF 0(3) 
model, the only new phase transitions generated were first order. We will 
consider in this letter the RP^ = S"^ / Z2 (real projective space) spin model 
in three dimensions. We place the spins on a cubic lattice with a nearest 
neighbors interaction: 

S = (5Y.{v,- v^)\ (1) 

<i3> 

where are normalized real three-components vectors. The local Z2 sym- 
metry Vi —>■ —Vi is preserved even after the SSB (Elitzur's theorem), and so, 
the sense of a spin is irrelevant, it is only its direction that matters. It is 
not hard to see 0, that in the AF case (|I|) is a lattice discretization of the 
action / tr[P(_R"^9^i?)^], where R E 5*0(3), and P is the diagonal matrix 
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{g,g, —g}, with g being the couphng. This is just a particular case of the 
NLcrM considered in reference 0. 

For j3 positive this model presents a weak first order phase transition 
which has been used to describe liquid crystals 0. The ordered phase cor- 
responds to states where all spins are aligned. 

For (3 negative there is also an AF ordered phase with a more complex 
structure p . There is a second order phase transition between the disordered 
phase and an ordered AF one. Let us call a site even (odd) when the sum 
of its coordinates x + y + z is even (odd). A state where, for instance, all 
spins on even sites are aligned in a given direction, and those on odd sites lie 
randomly in the orthogonal plane, has zero energy. So, at T = the ground 
state is highly degenerate with a global 0(2) symmetry. However, when 
fluctuations are taken into account, it can be shown in the low temperature 
limit , that both sublattices are aligned in mutually orthogonal directions, 
as a consequence of Villain's order from disorder mechanism . We remark 
that the remaining 0(2) symmetry is broken. We will show, with Monte 
Carlo simulations, that the breaking also holds in the critical region. 

In order to discuss the observables measured, let us construct the (trace- 
less) tensorial field T with components 

Tf = - (2) 

and its Fourier Transform TinaLxLxL lattice with periodic boundary 
conditions. 

The intensive staggered (non staggered) magnetization can be defined in 
terms of the tensorial field as the sum of the spins on even sites minus (plus) 
those on odd sites, or equivalently 

M, = :it(^,^,^) , (M = it(o,o,0)), (3) 

where V is the lattice volume. We have observed a phase transition at 
(3 ~ —2.41 for which is an order parameter (zero value in the disordered 
phase and a clear nonzero value in the L — > cxd limit in the ordered one). The 
magnetization M is also an order parameter. As these operators correspond 
to different irreducible representations of the translations group, we will study 
the scaling properties of each observable independently. 

To measure in a Monte Carlo simulation on a finite lattice we have con- 
structed scalars under the 0(3) group. For the magnetization and the sus- 
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ceptibility we compute respectively 



M = (VtrM2) , x = V (trM^) , (4) 



and analogously for the staggered observables. 

We have also measured the second momentum correlation length, which 
is expected to have the same scaling behavior at the critical point as the 



exponential (physical) one, but it is much easier to measure [jTO 

1/2 

Asin^{7i/L) 



where F is the mean value of the trace of T squared at minimal momentum 
{271 /L in direction x, y or z). To define we use Xs and compute Fg from 
T at momentum {2tt/L + vr, vr, vr) and permutations. 

The action (|I]) is suitable for cluster update methods by using the Wolff's 
embedding algorithm We have checked the performance of both the 



Swendsen-Wang and the Single Cluster methods [|Tl|]. Unfortunately, 
due to the AF character of the interaction, the critical slowing down is not 
reduced in any case, as there is always a large cluster that contains most of 
the lattice sites. In fact, we have measured a dynamic exponent z ^ 2 for 
both methods. 

We have also developed a Metropolis algorithm. Near the transition the 
spin fluctuations are large, and a spin proposal uniformly distributed over 
the sphere is accepted with nearly 30% probability. So we have used a 3 hits 
algorithm, reaching a mean 70% acceptance. 

Regarding the performance of the three methods mentioned for a given 
lattice size, the differences are very small in terms of the CPU time. We have 
selected the Metropolis method which is slightly faster. 

In table |l] we display the number of Monte Carlo sweeps performed for 
the different lattice sizes as well as the integrated autocorrelation times for 
the observables Xs and the energy. The total CPU time has been the 
equivalent of 12 months of DEC Alpha AXP3000 distributed over several 
workstations. 

Every 10 Monte Carlo sweeps we store individual measures of the en- 
ergy and of the Fourier transform of the tensorial magnetization at suitable 
momentum values. We have used the spectral density method [0 to ex- 



trapolate in a neighborhood of the critical point. The data presented here 
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22.14 


44.1 


21.4 


1.30 


24 


8.77 


107. 


48. 


1.82 


32 


10.13 


179. 


87. 


2.27 


48 


3.93 


430. 


205. 


3.10 



Table 1: Number of Monte Carlo sweeps performed for different lattice sizes. 
Measures have been taken every 10 sweeps. The integrated autocorrelation 
times (in sweeps) for both magnetizations and for the energy are also dis- 
played. The statistical errors are below the 5% level. We have discarded in 
each case about 200t^^ iterations for thermalization. 



correspond to simulations at two /? values (—2.41 and —2.4). We compute the 
quantities referred above as well as their /3-derivatives through the connected 
correlations with the energy. 

We have firstly analyzed several quantities that present a peak near the 
transition point. As we have found that the specific heat does not diverge, 
we have to limit ourselves to study quantities related with magnetization 
operators. The advantage of measuring a peak height is that its position also 
defines an apparent critical point allowing for a very simple and accurate 
measure. Unfortunately, quantities like the /9-derivatives of the magnetiza- 
tions or the connected susceptibilities (xcon = X~ VM'^) present their peaks 
far away from the critical point suffering from large corrections to scaling. 
For example, Xcon in the L = 16 lattice peaks aX (3 = —2.29 where the cor- 
relation length ^ is one half of its value at the critical point. This is due 
to the weakness of the tensorial ordering: the staggered magnetization, for 
instance, does not reach one half of its maximum until (5 < —3.5. 

Another possibility is to obtain the infinite volume critical point by other 
means and, then, to measure the different quantities at this point. By study- 
ing the matching of the Binder parameter for the staggered magnetization 
Vm, = 1- ((trM2)2)/(3(trM2)2), as well as that oi^s{L,/3)/L and the cor- 
responding non staggered quantities, we conclude that 

/?c e [-2.415,-2,405] . (6) 
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Figure 1: Deviations from a power law fit of the /3-derivative of the staggered 
(upper side) and non staggered (lower side) correlation lengths at the critical 
point, using data from lattice sizes with L > 16. The dotted lines correspond 
to a fit using all data sizes. 



To improve the above determination it is necessary a careful consideration 
of the corrections to scaling. This subject will be discussed in a forthcoming 
paper 0. 

In the case of quantities that change rapidly at the critical point as the 
magnetizations do, the errors in the determination of the critical point affect 
very much the results, and this method is not accurate. Nevertheless, we have 
found important quantities, like the /3-derivatives of the correlation lengths, 
which are very stable. Both d^s/dp and d^/djS should scale as L^^^f" . Fitting 
the data from all lattice sizes we obtain an acceptable fit: v = 0.793(2) with 
X^/dof = 2.0/5 and u = 0.787(2) with x^/dof = 5.5/5; however, if we discard 
the L = 6 data the fitted parameters change significantly and change again 
after discarding the L = 8 ones. The fits for L > 12 and L > 16 agree 
within errors. We, thus, choose the L > 12 data for computing u but take 
the statistical error from the fit with L > 16 (see figure |I|): 

: z/ = 0.788(7) 
e : z/ = 0.779(6). ^' 

To estimate the errors, not considered in (^, associated with the uncer- 
tainty in the determination of the critical point we repeat the fits with (3 at the 
limits of the interval @. We observe that u changes by an amount of a 1%. 

To avoid the problems reported above, we have also used a method di- 
rectly based on the Finite Size Scaling ansatz, that allows to write the mean 
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Table 2: Critical exponents obtained from a Finite Size Scaling analysis using 
data from lattices of sizes L and 2L. In the second row we show the operator 
used for each column. 



value of any operator O as 

(0(L,/3)) = LVo(e(A/?)/^) + ... , (8) 

where C,{L,P) is the correlation length measured at coupling /? in a size L 
lattice, fo is a smooth operator-dependent function and x depends also on 
O. The dots stand for corrections to scaling. 

Measuring (O) at the same coupling in lattices 2L and L and using (|[) 
we can write for their quotient 




Considering the dependence of Qo on p = ^(2L, /3)/^(L, /?), we just have 
to measure at the point where p = 2 to obtain Qo = 2^ up to corrections 
to scaling, x can be written in terms of the critical exponents: x = '-f/u 
for ^ = for ^) etc.. To compute z/ we can use O = d^/d(3 for 

which X = 1 + 1/v. As the quantity p is an observable and not an external 
parameter, like this procedure does not require a previous determination 
of I3c- Another advantage is that the result depends only on measures in just 
two lattices what allows a better error estimation. In columns 2 and 3 of 
table we report the results for the thermal exponent v obtained from data 
in lattices L (first column) and 2L, using as operators the /^-derivatives of 
the correlation lengths and ^ respectively. Even in the smaller lattices we 
do not observe any corrections to scaling. The data from columns 2 and 3 
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Figure 2: Quotients of several observables as a function of the quotient of 
the staggered correlation lengths. The sizes of the symbols are a growing 
function of the lattice size. is the /^-derivative of the staggered correlation 
length. 



are very correlated statistically and so, by taking the mean value, the errors 
are only slightly reduced. We select as our best estimation the mean of the 
results for the lattices 16-32: 

u = 0.784(8). (10) 

To be compared with the values z/(0(3)) = 0.704(6) and z/(0(4)) = 
0.748(9) [|1|]. 

In the case of the magnetic exponents, due to the large slope of Q as a 
function of p (see figure |^), a correct determination of the errors requires to 
take into account the statistical correlations of the whole data. We obtain for 
7s and jSg, in the most favorable cases, errors as small as a 0.1%. From these 
we compute, using the scaling relations 7s/z/ = 2 — r/^ and 2(3s/i' = D — 2 + ris, 



8 



the values of rjs with acceptable accuracy (see columns 4 and 5 of table 
In this case, the corrections to scaling are only significative for the L = 6 
lattice. We quote as our preferred value 

r]s = 0.038(2). (11) 

For the non staggered sector, we observe that the usual susceptibility 
diverges much more slowly than the staggered one [js — 1 ~ 1.02). The 
results for rj, using the corresponding scaling relations, are reported in the 
last two columns of table ^. In this case the corrections to scaling are non 
negligible for all lattice sizes. As the data fit very well to a linear function of 
1/L we take as the L ^ oo value 

7] = 1.338(8), (12) 

where the error is half statistical and half due to the possible deviations from 
linearity. Notice that the large value of r] means that, at the critical point, 
the spatial correlation function (G{r)\j3=i3^ ~ |r|^^^''), in the non staggered 
sector, decreases much faster than in the staggered case. 

Our results show that the 0(3) symmetry is fully broken in the ordered 
phase near the critical point. If we discard order terms {t being the reduced 
temperature), M is zero and the eigenvalues of M^, are {e, 0, -e} with e oc t^' 
{[3 — (3s ~ 0.51). Consequently, the magnetization tensors of the even and odd 
sublattices are opposite and the eigenvectors corresponding to the maximum 
eigenvalues are orthogonal, and similarly for the minimum ones. Considering 
the order terms the orthogonality will only hold approximately. 

We have studied a spin model in three dimensions with the symmetries 
of the 0(3) group but with very interesting new properties: it presents an 
ordered vacuum where the 0(3) symmetry is fully broken; the transition 
belongs to a new Universality Class, as the thermal exponent v is different 
from previously known and, finally, the model has two odd (magnetic type) 
relevant operators with different associated rj exponents. This may explain 
why previous perturbative calculations fail to work for this model. 

We think that in addition to the interest of the model by itself, the 
results suggest further studies of related models, like the addition of vector 
interactions, or four dimensional systems. 

We thank Alan Sokal for many enlightening discussions at the beginning 
of this work, specially regarding the structure of the vacuum in the T — > 
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Luis Alonso, Juan Jesiis Ruiz-Lorenzo and Alfonso Tarancon. This work 
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AEN93-0776. H. G. Ballesteros and V. Martin-Mayor are MEG fellows. 

References 

[1] T. Dombre and N. Read Phys. Rev. B39, 6797 (1989). 
[2] H. T. Diep, Phys. Rev. B39, 397 (1989). 

[3] P. Azaria, B. Delamotte and T. Jolicoeur, Phys. Rev. Lett. 64, 3175 
(1990); R Azaria, B. Delamotte, F. Delduc and T. Jolicoeur, Nucl. 
Phys. B408, 485 (1993). 

[4] M. Inui, S. Doniach and M. Gabay, Phys. Rev. B38, 6631 (1988). 

[5] J. L. Alonso, A. Tarancon, H. G. Ballesteros, L. A. Fernandez, 
V. Martin-Mayor and A. Muhoz Sudupe. To appear in Phys. Rev. B. 

[6] H. G. Ballesteros, L. A. Fernandez, V. Martin-Mayor and A. Muiioz 
Sudupe, work in preparation. 

[7] P. A. Lebwohl and G. Lasher, Phys. Rev. A6 426 (1972). 

[8] G. Korhing and R. E. Shrock, Nucl. Phys. B295 36 (1988); S. Romano, 
Int. J. of Mod. Phys. B8 3389 (1994). 

[9] J. Villain, R. Bidaux, J. P. Carton, and R. Conte, J. Piiys. (Paris) 41, 
1263 (1980). 

[10] S. Caracciolo, R. G. Edwards, A. Pehssetto and A. D. Sokal, Nucl Phys. 
B403, 475 (1993). 

[11] U. Wolff, Phys. Rev. Lett. 62, 3834 (1989). 

[12] R. H. Swendsen and J. S. Wang, Phys. Rev. Lett. 58, 86 (1987). 
[13] A. M. Ferrenberg and R. H. Swendsen, Phys. Rev. Lett. 61, 2635 (1988). 
[14] G. Holm and W. Janke, Phys. Lett. A173, 8 (1993). 
[15] K. Kanaya, S. Kaya, Phys.Rev. D51, 2404 (1995). 



10 



